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Abstract 



First measurements of the azimuthal anisotropy of neutral pions produced in PbPb 
collisions at a center-of-mass energy of a/snn = 2.76 TeV are presented. The am- 
plitudes of the second Fourier component (572) of the tt° azimuthal distributions are 
extracted using an event-plane technique. The values of are studied as a function 
of the neutral pion transverse momentum (pj) for different classes of collision central- 
ity in the kinematic range 1.6 < pj < 8.0GeV/c, within the pseudorapidity interval 
\rj\ < 0.8. The CMS measurements of ^(Pt) are similar to previously reported 71° 
azimuthal anisotropy results from a/snjv = 200 GeV AuAu collisions at RHIC, de- 
spite a factor of ~14 increase in the center-of-mass energy. In the momentum range 
2.5 < pj < 5.0 GeV/c, the neutral pion anisotropics are found to be smaller than those 
observed by CMS for inclusive charged particles. 
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A central goal of relativistic heavy-ion experiments is to create a deconfined phase of nuclear 
matter, the quark gluon plasma (QGP), at extreme temperatures and energy densities, and to 
characterize its properties. Observations at the Relativistic Heavy Ion Collider (RHIC) suggest 
that an extremely dense partonic medium with near-perfect fluid properties is formed HHII. 
These observations include the suppression of high-transverse-momentum (pj) hadron pro- 
duction, referred to as "jet-quenching"; strong azimuthal anisotropies in bulk particle produc- 
tion at low pj; and baryon-meson differences in hadron suppression patterns and azimuthal 
anisotropies at intermediate pj. Measurements of the azimuthal correlations of the produced 
particles play a key role in understanding the dominant physics processes in each of these 
transverse momentum ranges. 

At low pj (<2GeV/c), the azimuthal anisotropy of the emitted particles is understood to be 
the result of a collective hydrodynamic expansion of the medium, converting any initial-state 
spatial anisotropy (eccentricity of the nuclear overlap region) into a final-state momentum 
anisotropy ||5j|SJ. The strength of the anisotropy is characterized by the values of the Fourier 
coefficients, v n , of the expansion of the particle yields given by ^ cxli^ 2v n cos n((p — tp^p), 
where (p is the azimuthal angle of the outgoing particles and xpEF is the event plane angle recon- 
structed using the beam direction, z, and the azimuthal direction of the maximum transverse 
energy in each event. The second Fourier coefficient, Vi, is referred to as elliptic flow. At higher 
transverse momentum (pj > 6GeV/c), the azimuthal anisotropies have been attributed to the 
path-length dependence of energy loss in the medium due to the asymmetry in the reaction 
zone l|7lTPTT|. In the intermediate pj region, the RHIC data show an enhancement of baryon 
production ||T2llT3| and a larger V2 of baryons as compared to mesons |T4llT5| . This behavior 
has been interpreted as a signature of quark recombination as the dominant production mech- 
anism of moderate pj hadrons, which implies the existence of quark degrees of freedom in the 
medium produced at RHIC |T6|. 

The measurements of the elliptic anisotropy for inclusive charged particles, produced in PbPb 
collisions at a nucleon-nucleon center-of-mass energy of ^/snn = 2.76 TeV at the Large Hadron 
Collider (LHC), have been reported by ALICE H3, ATLAS EHL and CMS E0. The measured 
V2 coefficients exhibit similar strength and pj dependence as those measured at RHIC in AuAu 
collisions at y/smj = 200 GeV |TT7| . Differences in the elliptical flow of baryons and mesons 
can be indirectly tested via the comparison of the strength of the Vi signals for 7T° mesons and 
inclusive charged particles. 

This Letter presents the first measurement of elliptic flow of 7T° mesons as a function of pj in 
PbPb collisions at a center-of-mass energy of ^/s^n = 2.76 TeV. The data were recorded by 
the CMS experiment during the first LHC heavy-ion run in November 2010. The 7r° meson 
elliptic flow is measured in the pseudorapidity range \r\\ < 0.8, where tj is defined as rj = 
— In [tan (9/2)], and 6 is the polar angle between the particle momentum and the anticlockwise 
beam direction. The measurement is performed over the full azimuthal coverage < (p < 2n, 
and spans the range 1.6 < pj < 8.0GeV/c. 

The detectors used for this analysis are the barrel Electromagnetic Calorimeter (ECAL) and the 
Hadron Forward (HF) calorimeter, which have an rj acceptance of | tj \ < 1.4 and 2.9 < \t]\ < 5.2, 
respectively. Despite a wider pseudorapidity coverage of the barrel ECAL, these results are 
restricted to \rj\ < 0.8 in order to allow a direct comparison with the charged particle elliptic 
flow results [191. The barrel ECAL is located within a 3.8 T solenoidal magnetic field. The 
ECAL is made of lead-tungstate crystals that have a short radiation length (0.89 cm), and a 
small Moliere radius (2.19 cm). A more detailed description of the CMS experiment can be 
found elsewhere |20[ . 



2 



The minimum-bias event sample is collected using coincidences between the trigger signals 
from each side of the interaction point using the Beam Scintillation Counters (BSC) (3.23 < 
\t]\ < 4.65) or the HE Such a coincidence of minimum-bias trigger with bunches colliding in 
the interaction region suppresses any events due to noise, cosmic rays, out-of-time triggers, 
and beam backgrounds. The trigger accepts (97 ± 3)% of the total inelastic PbPb cross section. 
Collision centrality defined as the fraction of total inelastic nucleus-nucleus cross section, is 
calculated using the sum of transverse energy (Ej) in towers from HF at both positive and 
negative z positions [21J. In this Letter, we present results based on centrality intervals of 10% 
width, ranging from 20-30% (more central) to 70-80% (more peripheral). For the most central 
collisions (0-20%), a small signal-to-background ratio limits the identification of tt° mesons. 

The 7T° mesons are measured by reconstructing their decay photons (tz° — > 77) in the bar- 
rel ECAL. Electromagnetic showers are found in the ECAL by forming clusters of contiguous 
crystals with a seed crystal having energy above a threshold of 200 MeV. Clusters are identi- 
fied as photons on the basis of a shower shape requirement called the S4/S9 ratio. Photons 
are reconstructed using a 3 x 3 array of crystals, which contain on average 93% of the photon 
energy. The quantity S4 is the total energy in a 2 x 2 array (a sub-matrix of the 3x3 array) 
containing the crystal with the highest energy deposited and S9 is the total energy in the 3x3 
crystal matrix. There are four possible 2x2 matrix combinations, and S4 is defined as the 
most energetic of these four combinations. Clusters with S9 > 400 MeV and S4/S9 > 0.87 are 
selected as photon candidates for tt° meson reconstructed invariant mass, m 7l7j calculations. 
The invariant mass of a photon pair (7,, jj) as measured in the ECAL is calculated from the 



energies and positions of the clusters, as given by m li7j = w2E,-£y(l — cos fly), where is the 

opening angle between the two photons. Candidate pairs are formed from each photon cluster 
in an event in a particular pj bin for the n° meson invariant mass calculation. A pj-dependent 
opening angle requirement and a cluster pair separation cut are also applied to the tc° meson 
invariant mass distribution. Pairs are selected with > -j^ + where a and b are opening 

angle cut parameters obtained from a detailed PYTHIA 6.422 simulation l|22| . The values of the 
parameters a and b are 0.17GeV/c and —0.11 (GeV/c) 2 , respectively. Further, a photon pair is 
rejected if the separation of the two photon clusters (distance is calculated based on the rj and 
(p coordinates of the clusters and using 1.29 m radius of the ECAL) is less than a threshold dis- 
tance at a certain pj. The threshold distance between photon-clusters decreases monotonically 
from 15 cm at pj « 1.6 GeV/c to 5.0 cm at pi ~ 8.0 GeV/c. At sufficiently high pj, photons from 
a nearly symmetric decay (E 7l ~ E 72 , where E 7i is the energy of a photon) can produce showers 
in the calorimeter that are reconstructed as a single cluster. In CMS, this effect is first visible 
around pj > 8.0GeV/c. Consequently, results presented here are restricted to pj < 8.0GeV/c. 

The 7T° meson yields are extracted statistically by subtracting the combinatorial background 
from the rc° candidate invariant mass distribution. The combinatorial background is estimated 
and subtracted using an event mixing technique, which forms pairs from photon candidates in 
different events. Each photon candidate is combined with all other photon candidates in three 
other events. The mixing of events is performed within intervals of centrality, z-vertex posi- 
tion and the event-plane angle fl9| orientation to replicate the background from uncorrelated 
pairs. All selections applied to the combinations of same-event pairs are also applied to mixed- 
event pairs. Event mixing is done in six z-vertex intervals of width Az = 5.0 cm in the range 
|z| < 15 cm. Similarly, the event-plane angle [19] is also divided into six intervals in the range 
< ipEF < Tt. The event-plane angle is determined from the HF, with flattening and resolution 
correction factors applied as in |T9| . The 7i° reconstruction efficiencies as a function of pj, cen- 
trality, and event-plane are studied by embedding simulated 7T° mesons in real events. A total 
of 100 k such events are analyzed, where each event has ten zr mesons embedded with a flat 
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pj and (p distribution over a range of 0.2 < pj < 10.0 GeV/c and \tj\ < 1.0 to avoid any edge ef- 
fects. The results for 7T° meson elliptic flow are corrected for the dependence of the reconstruc- 
tion efficiency on pj for all centralities. In addition an in-plane versus out-of-plane dependence 
is observed for the tc° meson reconstruction efficiency in more central collisions. Corrections 
for this effect range from 16% (1.6 < p T < 2.0GeV/c) to 6% (2.5 < p T < 3.0GeV/c) for the 
20-30% centrality interval. For higher pj intervals in this centrality class, such (^-dependent 
efficiency corrections are not needed. Similarly for the 40-50% and 50-60% centrality intervals, 
these corrections range from 7% (1.6 < p T < 2.0GeV/c) to 4% (2.5 < p T < 3.0GeV/c), while no 
^-dependent efficiency corrections are needed for the more peripheral events. 

Figure [l] (top panel) presents the 71° meson invariant mass distribution before background sub- 
traction for 2.5 < pj < 3.0GeV/c for the 40-50% centrality interval. This panel shows the 
same-event distribution (solid circles) and the mixed-event normalized background (dashed 
line). For a given pj bin, the mixed-event background distribution is normalized to the same- 
event signal distribution in the range 200-250 MeV/c. Different normalization regions such as 
175-225 MeV/c and 225-275 MeV/c are also studied and no significant change in the resulting 
7T° meson is observed. The middle panel of Fig. [l] shows the combinatorial-background- 
subtracted 71° meson invariant mass distribution (solid circles) in the same pj bin and central- 
ity class. Over-subtraction is observed for higher mass regions, m 77 > 250MeV/c. Investiga- 
tions using pythia and hydjet (1.8) |23| simulations show that this effect can be attributed 
to a correlated conversion background (converted photons) which has a different shape than 
a purely combinatorial background. By definition, the event-mixing technique cannot account 
for the effect of a correlated conversion background. Open symbols in the middle panel cor- 
respond to HYDJET simulations, the result obtained without rejecting any converted photons. 
The background-subtracted mass spectrum predicted by simulations is seen to reproduce the 
data well. HYDJET simulation results also show that the over-subtraction at high invariant mass 
is eliminated when the clusters from the converted photons are suppressed, as shown in the 
bottom panel of Fig. [I] The event yield is calculated by integrating the data in a two standard 
deviations (a, in units of mass) window around the mean (p) of the distribution. The a and 
p are determined from a Gaussian fit to the combinatorial-background-subtracted zr meson 
invariant mass distribution for every pj and centrality interval. To avoid any model depen- 
dence, no corrections to the data are applied in order to account for these converted photons; 
instead asymmetric mass integration ranges of p — 2a < m 77 < p, and p — 3a < m 77 < p are 
employed to understand the systematic effect of the conversion contribution to the mass peak 
in the higher mass regions. Studies showed that the maximum effect of the correlated back- 
ground on the yield extraction in the mass integration range is less than 16%. To obtain the 
dependence of n° meson production on azimuthal angle, the extracted yield is first measured 
in a given pj bin as a function of the azimuthal angle between the n° meson trajectory and the 
event-plane orientation, \pYp, found as described in Ref. [19 j. The measurement is performed 
in six equally spaced intervals of Acp = <P(tz°) — t/?EP in the range < Acp < n/2. The 7T° meson 
yields corrected for reconstruction efficiency are measured for each Acp bin and the resulting 
angular distribution, dN / dA<p, is fitted with No (1 + 2v2 cos 2A<p) to determine the strength of 
the modulation in the yield. We use an analytic linear % 2 fitting procedure that matches the 
integral of No (1 + 2??2 cos 2A</>) over each A(p bin to the measured n° meson yield within the 
corresponding bin |Tl4l[T5| . 

Systematic uncertainties are assessed by varying the S4/S9 ratio and the mass integration 
ranges. A combination of the S4/S9 = 0.87 and |m 77 — p\ < 2.0cr mass integration range serves 
as a reference in this analysis. The n° meson V2 results are calculated for S4/S9 = 0.83 or 
0.91 keeping the mass integration range at a reference value of |wj 77 — p\ < 2.0a. In addi- 
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Figure 1: (Color online) Top panel: 7T° meson invariant mass distribution with combinato- 
rial background for 2.5 < pj < 3.0GeV/c for 40-50% centrality interval in PbPb collisions 
at a/snn = 2.76 TeV. The solid line represents the combinatorial background estimated from 
the event-mixing technique. Middle panel: Combinatorial background-subtracted n° invariant 
mass distribution (solid circles) in the same pj bin and centrality interval. Open squares show 
HYDJET simulation results. Bottom panel: HYDJET simulation results obtained after rejecting 
converted photons. The HYDJET simulation results have been scaled up by a factor of 556 to 
match the data. 

tion to asymmetric mass integration ranges, symmetric ranges such as |m 77 — \i\ < 3.0a, and 
l m 7? ~ Y-\ < 1.5(7" are used to determine the n° meson vi results for all centralities keeping 
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S4/S9 fixed at 0.87. The largest observed differences in the V2 results based on different S4/S9 
ratio cuts and m 77 — p ranges are used to determine the systematic uncertainty. The systematic 
uncertainty determined from the precision of the ^-efficiency curves obtained from the embed- 
ding procedure ranges from 18% to 4% from the lowest to the highest px intervals for 20-30% 
centrality. For 70-80% centrality the systematic uncertainty varies from 7.2% to 9%. The total 
systematic uncertainties obtained upon adding all the sources listed above in quadrature vary 
from 21% (1.6 < p T < 2.0 GeV/c) to 31% (6.0 < p T < 8.0 GeV/c) for the 20-30% centrality inter- 
val. Similarly for 70-80% these uncertainties change from 9.6% (1.6 < p T < 2.0GeV/c) to 33% 
(6.0 < px < 8.0GeV/c). Systematic uncertainties arising from the trigger efficiency are found 
to be negligible. 

The 7i° meson »2(Pt) results are shown in Fig. |2]for six centrality classes from 20-30% to 70- 
80%. CMS 7T° meson V2 results, shown as solid circles, are compared to PHENIX n° V2 results 
0, for AuAu collisions at -^s^iv = 200 GeV, shown as open circles. Our measurement shows 
qualitatively similar features as observed at RHIC energies despite an order of magnitude in- 
crease in the center-of-mass energy. This observation is consistent with elliptic flow results for 
charged particles at RHIC and LHC fTTllTBll . Green (grey) shaded bands show the systematic 
uncertainties associated with the CMS n meson (charged particle) measurements. 

Figure [2]also presents a comparison between CMS 7r° meson V2 results (solid circles), and CMS 
inclusive charged particle vi [19[ (open squares) as a function of px using the event-plane 
method. The 7r° meson V2 is systematically lower than that for inclusive charged particles 
V2 between 2.5 < px < 5.0GeV/c for all six collision centrality intervals. For particles with 
intermediate px at RHIC, the vi values of baryons are observed to be higher than those for 
mesons |l"4"llT5| . The differences observed between the inclusive charged particle and 7r° me- 
son results may be due to the contribution from baryons which would increase the overall vi of 
the inclusive charged particles, compared to that for neutral pions, assuming a baryon-meson 
V2 splitting comparable to that seen at RHIC. 

In summary, the CMS detector has been used to perform the first measurements of the az- 
imuthal anisotropy of neutral pions in PbPb collisions at y/s^ = 2.76 TeV. The measure- 
ments of V2 were presented as a function of px for six centralities, from 20-30% to 70-80% for 
1.6 < px < 8.0GeV/c. Results were compared with PHENIX tt° meson [9| and CMS inclusive 
charged particle measurements |Tl9l . It was found that the values of z?2(Pt) for neutral pions 
measured at RHIC and the LHC were of comparable magnitude. In the momentum range 
2.5 < px < 5.0 GeV/c, the magnitude of elliptic flow for neutral pions was found to be system- 
atically lower than that for charged particles. This behavior is consistent with observations at 
lower collision energies, where this difference is found to be caused by the larger elliptic flow 
of baryons compared to mesons. 
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Figure 2: (Color online) CMS tt° meson Vi (solid circles) compared to PHENIX 7T° meson (9l 
(open circles) for mid-rapidity ( | rj \ < 0.8 and \tj\ < 0.35, respectively) and CMS charged particle 
V2 (open squares, \rj \ < 0.8). Results are presented as a function of pj for six centrality intervals 
(20-30% to 70-80%). Green (grey) shaded bands represent systematic uncertainties associated 
with CMS 7i° meson (charged particle) measurements. Only statistical uncertainties are 
shown for the PHENIX results. The systematic uncertainties for the 50-60% centrality on the 
PHENIX data points are 10.4%. 
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